Gate control of a quantum dot single-electron spin in realistic confining potentials: 

anisotropy effects 
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Among recent proposals for next-generation, non-charge-based logic is the notion that a single 
electron can be trapped and its spin can be manipulated through the application of gate potentials. 
In this paper, we present numerical simulations of such spins in single electron devices for realistic 
(asymmetric) confining potentials in two-dimensional electrostatically confined quantum dots. Using 
analytical and numerical techniques we show that breaking the in-plane rotational symmetry of the 
confining potential leads to a significant effect on the tunability of the g-factor with applied gate 
potentials. In particular, anisotropy extends the range of tunability to larger quantum dots. 
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I. INTRODUCTION 

The notion of using single-electron spins for quan- 
tum computing and next-generation logic is an attractive 
idea that has received considerable attention in recent 
years . 1 - 2 - 3 - 4 - 5 - 6 - 7 - 8 - 9 - 10 - 11 - 12 - 13 

In order to integrate new concepts with existing 
semiconductor technology, a number of researchers 
have recently explored the possibility of using elec- 
tric fields generated by imposed gate potentials to 
manipulate single-electron spins in quantum dot de- 



vices. 



14, 15, 16, 17, 18, 19, 20. 21, 22. 23, 24, 25. 26. 27. 28. 29, 30, 31, 32. 3:5. 34 



The goal of the present work is to utilize state-of- 
the-art numerical techniques to explore the fundamental 
physics of single-electron spin devices and to provide real- 
istic information for the practical design of such systems. 
We utilize a finite-element based numerical technique to 
study electrostatically defined quantum dots that is sim- 
ilar to other recently published work. 35 - 

A key result of the present work is the discovery 
that spatial symmetry breakin g 36 ' 37 i 38 resulting from the 
anisotropy of realistic confining potentials results in a 
significant enhancement of the electric-field tunability of 
the electron <?-factor over that found for symmetric po- 
tentials. 



II. COMPUTATIONAL METHOD 

We utilize a multi-scale, multi-physics simulation strat- 
egy based on the finite element method 3 ^ to provide a re- 
alistic description of the physics of single-spin devices in 
three-dimensional geometries. The ideal is to solve self- 
consistently the Maxwell equations of electrostatics with 
the Schroedinger equation in three dimensional geometry. 
Unfortunately such a solution is not feasible given cur- 
rently available techniques due to the disparity of length 
scales in the problem. We thus seek an approximate so- 
lution that is built up in stages. 

In the first step of our approach, we construct a three- 
dimensional model of the device and calculate the gate- 
induced electrostatic potentials that cause the formation 
of a quantum dot in the two-dimensional electron gas 
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FIG. 1: (color online) Electrostatic potential for a prototype 
single-electron device plotted in the 2DEG layer. This fig- 
ure illustrates a single-spin device consisting of two triangular 
gates above a 2DEG. The gates were held at IV and the 2DEG 
was held at OV. For simplicity of the electrostatic calculation, 
the 2DEG was treated as a classical perfect conductor. The 
dimensions of the device in the x and y directions are 2.8/im 
and 1.8/xm respectively and the thickness is 



(2DEG) at a AlGaAs/GaAs heterojunction as illustrated 
in Figs.Q] [21 and[3] This geometry corresponds to proto- 
type devices that are under consideration by experimen- 
talists at the University at Albany, State University of 
New York. 

In order to obtain the electrostatic solution for the con- 
fining potential, we approximate the 2DEG as a classical 
perfect conductor and give it a finite width. The width 
(sa 0.05/im) is unrealistically large from a quantum per- 
spective but is assumed to give a reasonable description 
of the spatial variation of the potential in the layer of the 
2DEG. In a subsequent step we treat the 2DEG from a 
realistic quantum mechanical perspective^ 

Figs. [2] and [3] are one-dimensional plots obtained from 
Fig. [T] by plotting the potential along a line in the 2DEG 
along high symmetry directions. These one-dimensional 
potentials are then fit to polynomial forms P x {x) and 
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FIG. 2: Electrostatic confining potential in the 2DEG along 
the symmetry axis of a prototype single-electron device. This 
figure was made by plotting the potential of Fig. [T] along a 
line in the 2DEG through the symmetry axis of the device 
(the x-axis of Fig. [1] i.e. a line running from one gate to the 
other) intersecting with the central region. 
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FIG. 4: Heterojunction self-consistent potential and lowest 
two wave functions of the 2DEG. These results demonstrate 
consistency of our results with other published work.— 
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FIG. 3: Electrostatic confining potential in the 2DEG normal 
to the symmetry axis of a prototype single-electron device. 
This figure was made by plotting the potential of Fig. [T] along 
a line in the 2DEG normal to the symmetry axis of the device 
(the y-axis of Fig. [1} and intersecting the central region. 

P y (y) . These are then used as a potential of the form 

V real (x,y) n P x (x) + P y (y) (1) 

to approximate the confining potential of the electron 
in the Schroedinger equation. Before considering elec- 
tron motion in the above potential, V rea i, we consider 
the simpler quadratic potential 

V qU ad = ^mcu 2 (ax 2 + f3y 2 ) (2) 

that allows for systematic studies. For convenience we 
have written the strength of the potential in harmonic 
oscillator form by defining the pre- factor ^muj 2 . 



In the second step we calculate the wave functions and 
self-consistent potential at the heterojunction between 
AlGaAs and GaAs that describes the formation of the 
two dimensional electron gas as illustrated in Fig. [4] We 
do this calculation primarily to benchmark our numerical 
method by making contact with a well-know result from 
the literature^ 

The results of Fig. 2] are obtained by solving 
the following coupled equations that constitute the 
self-consistent Schroedinger-Poisson equations including 
exchange-correlation effects 

-h 2 d ( 1 dibJz)\ , N , N 

T£UW 4 + v{z)Mz) = EiMz) (3) 
Tz{ £oK{z) ^f) = 'Y.OiM')? - fto (4) 

where k[z) and p{z) are the fixed spatially dependent 
dielectric function and background charge density of the 
the interface as described in Ref. |40|. The potential en- 
ergy of the 2DEG is given by: 

V(z) = -e<t>{z) + V xc {z) (5) 

where 4>{z) is the self-consistent potential and V xc (z) is 
the exchange-correlation potential. In the above equa- 
tions, the coordinate z is measured relative to the inter- 
face between AlGaAs and GaAs. 

The results of Fig. [4] are in excellent agreement with 
previous results 40 confirming the soundness of our ap- 
proach. 

We next consider the formation of an electrostatically 
defined quantum dot by applying a symmetric, confining 
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FIG. 5: (color online) Quantum dot wave function plotted 
in the x-z plane formed by applying a quadratic confining 
potential in the plane (i.e. along the x axis). This potential 
is characterized by the parameter £q — 20nm (quantum dot 
radius, see Eg. I16|l . 



potential in the plane of the 2DEG as illustrated in Fig. [5] 
In other words, we add a potential of the form 



V x (x) 



1 



2 2 

mu) n x 



(0) 



to Eq. [5] and solve the system of Eqs. [U 0] and [5] self- 
consistently in the two-dimensional x — z domain. Fig- 
ure [5] clearly illustrates the formation of a quantum dot 
in the potential well of the 2DEG as expected. 

In the remainder of this paper we focus our attention 
on motion in the plane of the 2DEG and contrast ef- 
fects associated with quantum dots in symmetric and 
asymmetric confining potentials as illustrated in Figs. [6] 
andEl respectively. These figures were obtained using the 
quadratic model potential of Eq. [5] with a — (3 — 1 and 
a = 1; f3 = 2.8 respectively in the two-dimensional in- 
plane (i.e. x— y plane) Hamiltonian, H xy , to be discussed 
in the following. The parameters of the asymmetric po- 
tential were chosen so as to mimic the realistic potential 
of Figs. [2] and [3l The wave function in the asymmetric 
model potential of Fig. [7] should be contrasted with the 
wave function in the realistic potential (i.e. using the 
form of Eq. [TJ) as shown in Fig. [5] 

We consider the motion of the electron in the x — y 
plane of the quantum dot in the presence of a magnetic 
field oriented perpendicular to the plane of the 2DEG. 
Our approach closely follows that of Ref. [IJ. Thus the 
total Hamiltonian can be written as 



H — H xy + H z + H s 



(7) 



where H z corresponds to motion normal to the interface 
(as discussed in the context of Eqs. [21 [HandO, H so is 
the spin-orbit interaction to be discussed shortly and the 
remaining term is given by: 
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+ -muj (ax + fiy ) + g jjL B o z B, (8) 



FIG. 6: (color online) In-plane wave-function for quantum dot 
formed by a symmetric quadratic confining potential (Eq. [2} 
in the (x — y) plane with £ a = 40nm (see Eq. I16p . 
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FIG. 7: (color online) In-plane wave-function for quantum dot 
formed by an asymmetric quadratic confining potential in the 
(x — y) plane. The parameters for this potential where chosen 
to mimic the realistic potential of Figs. [2] and [3] respectively 
and are given by l = 30nm (see Eq. I16p and a = 1 and 
P = 2.8 (see Eq.©. 



where the kinetic momentum operator: 



P = p+-A, 

c 



(9) 



is the sum of the canonical momentum 

p = -ih{d x ,d y ,o), (10) 

and the vector potential (in the symmetric gauge) 



(ii) 



The eigenstates of H xy (Eq. [5]) with a — (3 are the well- 
known Fock-Darwin states i 41 ' 42 The situation with q//5 
also has an analytic solution^ We have verified that our 
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FIG. 8: (color online) In-plane wave-function for quantum dot 
formed by the potential of Eq. [1] illustrated in Figs. [2] and [3] 



numerical solution of H xy \ip >= e|-0 > is consistent with 
these analytical results. 

Lastly we consider the spin orbit interaction as em- 
bodied in the Hamiltonian H so which is the essential in- 
gredient in the phenomena of electric field induced spin 
switching We write: 

H so = H R + H m + H D2 (12) 

where the Rashba interactio n 44 ' 45 is given by: 



Hr = (vxPy - (TyPx \ , 



(13) 



and the linear and cubic Dresselhous interaction a 46 ' 47 are 
written as: 

H m = ° J ™ 7cfc2 (-a x P x + a y P y ^j , (14) 

which is linear in components of the momentum operator 
P and 



TT _1C 

tlDI — To- 



-a x P x P' 



VyPyP* 



(15) 



which is cubic in components of the momentum operator 
(h.c. denotes the Hermitian conjugate). 14 Note, the elec- 
tric field strength E that enters Eq. [T3] is that associated 
with the heterojunction \E\ = dV{z)/dz and is treated as 
an adjustable parameter. Physically we can implement 
changes in E through the application of appropriate gate 
potentials. All numerical parameters in the above pieces 
of H so are those for GaAs found in Ref. [bj 

The eigenvalue equation H \ tp > = e \ ip >, with H given 
by Eqs. [7] through [151 was solved numerically to obtain 
the lowest few eigenstates and eigenenergies vs. the vari- 
ous parameters of the system. These parameters include 
the magnetic field strength B, the electric field E, and 



the strength of the quantum dot confinement potential. 
This latter parameter is conveniently characterized by 
defining the quantum dot radius as follows: 



mw 



(16) 



The notion of electric field induced spin switching is 
quantified by defining an effective electron g factor by 
the following definition: 



(17) 



to describe the energy difference between the lowest en- 
ergy up and down spin states. Thus we consider the 
lowest two states (including spin) e 2 and t\ and calculate 
the effective g factor as: 



(£2 -d) 
[i B P 



(18) 



Results for the variation of this effective g factor as a 
function of the parameters E, B and t D is presented in 
the following section. 



III. RESULTS 

We now turn to a presentation of the key results of this 
work: the tunability of the electron (/-factor through the 
application of electric and magnetic fields. 

Figure [5] is consistent with previous published worki 4 ^ 
and illustrates the g-factor tunability vs. the strength 
of the applied electric field and confining potential (as 
parametrized by the quantum dot radius £ a ) for fixed 
magnetic field (B = IT) for the symmetric quantum dot 
in the quadratic potential of Eq. [2] with a — j3 = 1 . 

Figure [10] is also consistent with previous published 
workii and illustrates the g-factor tunability vs. the 
strength of the applied electric field and and magnetic 
field for fixed confining potential (parametrized by the 
quantum dot radius £ a — 20nm). 

Upon introducing in-plane anisotropy to the confining 
potential we find significant changes in the electric-field 
induced g-factor tunability as illustrated in Fig.[TTJ This 
figure was generated by choosing a = 1, /3 = 2 and l Q = 
120nm. 

To quantify the effects of in-plane anisotropy, we have 
carried out a parameter study of the g-factor tunability 
vs. degree of anisotropy and the results are presented in 
Fig. 1121 This figure was generated by fixing the quan- 
tum dot radius at £ a = 120nm and holding a = 1 while 
varying (3 with B = IT. 

Lastly we consider the results for g-factor tunability for 
quantum dots in the realistic potential of Figs. [5] and [3] 
and Eq.[TJ Figure [T"3l illustrates the results for a quantum 
dot in the realistic potential (middle curve) in compar- 
ison with symmetric (lower curve) and asymmetric (up- 
per curve) quantum dots. The parameters were chosen 
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FIG. 9: (color online) Electric field induced changes in the 
g-factor vs. quantum dot radius for various electric field 
strengths for a symmetric quantum dot in the quadratic po- 
tential of Eq. [2] From top to bottom, the curves represent in- 
creasing electric field strength as follows. The first curve cor- 
responds to 1 x 10 4 V/cm, and the rest range from 1 x 10 5 V/cm 
through 1 x 10 6 V/cm in equal steps with B — IT. This results 
is consistent with Ref. 14. 
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FIG. 10: (color online) Electric field induced changes in the g- 
factor vs. magnetic field for various electric field strengths for 
the symmetric quantum dot. From top to bottom, the curves 
represent increasing electric field strength as follows. The first 
curve corresponds to 1 x W 4 V/cm, and the rest range from 
1 x W 5 V/cm through 1 x l0 6 V/cm in equal steps. For this 
calculation, the quantum dot radius was fixed at £ = 20nm 



to closely mimic the realistic potential and are given by 
£ = 30nm and a = (3 — 1 for the symmetric potential 
and a = 1; /3 = 2.8 for the asymmetric potential. 

From Fig. [13] we conclude that the realistic potential 
result is bracketed by those for symmetric and asymmet- 
ric model potentials and, as such, contains characteristics 
of each. 
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FIG. 11: (color online) Electric field induced changes in the 
g-factor vs. quantum dot radius for various electric field 
strengths for an asymmetric quantum dot (wave function not 
shown). From top to bottom, the curves represent increasing 
electric field strength as follows. The first curve corresponds 
to 1 x 10 4 V/cm, and the rest range from 1 x 10 5 V/ cm through 
1 x 10 6 V/ cm in equal steps. Again we choose B — IT. 
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FIG. 12: (color online) Electic field induced changes in the g- 
factor vs. the degree of anisotropy of the quantum dot confine- 
ment potential for various electric field strengths. From top to 
bottom, the curves represent increasing electric field strength 
as follows. The first curve corresponds to 1 x 10 4 V/ cm, and 
the rest range from 1 x 10 5 V/cm through 1 x 10 6 V/cm in 
equal steps. 



IV. DISCUSSION 

The key result of this work is illustrated in Figs. [5] 
and llll anisotropy in the confining potential significantly 
extends the size range of quantum dots that exhibit elec- 
tric field induced <?-factor tunability. Indeed, in Fig. [9] 
we see that all of the curves collapse onto a single curve 
for large quantum dots (i.e. larger than £ a — 120nm) 
negating the switching effect. With anisotropy, however, 
there is no degradation of the switching effect for large 
dots. 

Another important result of this work is the realization 
that the degree of anisotropy need not be very large in 
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FIG. 13: (color online) Changes in the g-factor vs. electric 
field for quantum dots in the realistic potential of Figs. [2] 
and (using the form of Eq. [TJ in comparison to results for 
symmetric and asymmetric model potentials. The black (mid- 
dle) curve corresponds to the realistic potential while the red 
(lower) and green (upper) curves correspond to symmetric 
and asymmetric model potentials respectively. The parame- 
ters were chosen to closely mimic the realistic potential and 
are given by l = 30nm and a = f3 = 1 for the symmetric 
potential and a = 1; j3 = 2.8 for the asymmetric potential. 

order to obtain significant changes in the gate induced 
g-factor tunability, as illustrated in Fig. [T3J We see that 
the maximum effect is obtained for roughly f3/a ~ 1.3 
and begins to decrease slightly for larger values of this 
shape-anisotropy ratio. The jumps in the value of the 
g-factor from positive to negative value are indicative of 
level crossings (e.g. the relative ordering of spin up and 



down levels changes as a function of the anisotropy). 

Lastly we have seen from Fig. Q2] that results for quan- 
tum dots in realistic potentials contain characteristics of 
both the symmetric and asymmetric model potentials. 
This emphasizes the need for realistic numerical simula- 
tions since the physics of such systems cannot be fully 
captured by symmetric or asymmetric model potentials 
alone. 



V. CONCLUSIONS 

We have carried out a numerical simulation study of 
gate induced tunability of the electron g-factor in a pro- 
totype single electron spintronic device. We consider a 
realistic three-dimensional geometry and employ a nu- 
merical approach based on the finite element method 39 . 

Due to the large disparity in physical length scales we 
have adopted an approximate strategy as a complete self- 
consistent solution of the full problem is prohibitive given 
existing computational tools. In our approach we have 
investigated the problem in stages. 

We find that symmetry breaking due to anisotropy of 
the quantum dot confining potential leads to significant 
changes in the gate induced tunability of the g-factor vs. 
quantum dot radius extending the size range of quantum 
dots that can be tuned. We also find that the anisotropy 
ratio need only be on the order of /3/cv = 1.3 to produce 
the largest effects. 

By employing non-perturbative, fully numerical meth- 
ods and realistic geometries, our approach is providing 
insights that might be difficult or impossible to obtain 
using analytical techniques alone. 
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